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The Green Propellant Infusion Mission (GPIM) will demonstrate the capability of a green 
propulsion system, specifically, one using the monopropellant, AF-M315E. One of the risks identified 
for GPIM is potential contamination of sensitive areas of the spacecraft from the effluents in the 
plumes of AF-M315E thrusters. Plume characterization of a laboratory -model 22 N thruster via optical 
diagnostics was conducted at NASA GRC in a space -simulated environment. A high-frequency pulsed 
laser was coupled with an electron-multiplied ICCD camera to perform Raman spectroscopy in the 
near-field, low-pressure plume. The Raman data yielded plume constituents and temperatures over a 
range of thruster chamber pressures and as a function of thruster (catalyst) operating time. Schlieren 
images of the near-field plume enabled calculation of plume velocities and revealed general plume 
structure of the otherwise invisible plume. The measured velocities are compared to those predicted 
by a two-dimensional, kinetic model. Trends in data and numerical results are presented from catalyst 
mid-life to end-of-life. The results of this investigation were coupled with the Raman and Schlieren 
data to provide an anchor for plume impingement analysis presented in a companion paper. The 
results of both analyses will be used to improve understanding of the nature of AF -M315E plumes and 
their impacts to GPIM and other future missions. 

I. Introduction 

T HE Green Propellant Infusion Mission (GPIM) is a Technology Demonstration Mission (TDM) project, 
sponsored by NASA’s Space Technology Mission Directorate (STMD). 1 The goal of GPIM is to 
demonstrate the capability of a “green” propulsion system, specifically, one using the monopropellant 
AF-M315E. Green monopropellants offer lower toxicity and reduced safety hazards, while also aiming to 
provide higher performance than the current state of the art fuels. The GPIM propulsion system will be flown 
as a payload on a Ball Aerospace BCP-100, a small, standardized spacecraft. The propulsion system was to 
have one 22 N thruster for primary divert maneuvers and four 1 N thrusters for attitude control. 
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One of the risks identified for GPIM is potential contamination of sensitive areas of the spacecraft from 
the effluents in the plumes of AF-M315E thrusters. The plume risk mitigation activities include modeling 
the plume flow fields of the AF-M315E thrusters, assessing the plume impingement on the BCP-100 
spacecraft including the impact on the power generating capabilities of the solar arrays, and conducting 
ground-based plume measurements on an AF-M315E thruster to correlate the plume modeling with plume 
data. 

The Raman diagnostic was designed to interrogate the plume of a high-expansion area nozzle with exit 
plane pressures on the order of 0. 1 psia. This environment has proven very difficult to achieve Raman signals 
in the past because the signal strengths are proportional to the density. 2,3 Most Raman combustion diagnostics 
are implemented at or above atmospheric pressure suggesting a two or three order of magnitude decrease in 
signal using conventional approaches. In addition to significantly lowering the pressure, nozzle expansion 
also reduces the temperature of the target gas well below where Raman diagnostics typically interrogate 
combustion. This removes many of the ro -vibrational characteristics that are used to quantify both 
temperature and species concentrations. In previous experiments, high-power pulsed lasers were required to 
achieve modest signal to noise ratios and focusing these lasers to a very small interrogation volume posed a 
significant hazard to the windows of the vacuum chamber. In this investigation, a high-frequency laser was 
coupled with an electron-multiplied ICCD (EMCCD) camera to mitigate the constraints imposed by high- 
power, single-shot lasers and allows accumulation of multiple shots on the ICCD in preference to ensemble 
averaging. 

The objectives of this investigation include identifying the major constituents in the AF-M315E plume 
and correlating them to variations in chamber pressure, time on catalyst, and spatial position in the plume. 
These data will be included in a plume-impingement analysis. 4 Also required for that analysis are velocities 
in the plume which will obtained via schlieren imaging. The images may also reveal any significant radial 
variation in plume structure associated with non-uniform combustion. This investigation also serves as the 
demonstration of capability for the diagnostics needed to interrogate plumes from higher, more flight -like, 
expansion ratio nozzles with AF-M315E or other green propellants. 

II. Theoretical Background 


The following overview of Raman scattering theory, Schlieren diagnostic technique and VIPER 
modeling is intended to facilitate understanding of the data and calculations discussed in the next section. 
More thorough presentations are available in the literature and are noted where applicable. 

A. Raman Diagnostic 


The theory of Raman scattering in diatomic and polyatomic molecules is well established. 5 According 
to Placzek’s polarization theory of the Raman effect, the corresponding differential Raman scattering cross- 
section per molecule at thermal equilibrium is given by 5 
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Here, the superscript 6 indicates the relative angle of the planes of polarization between scattered and 
incident light; the subscripts of i and f stand for the initial and final states respectively; s 0 , is the vacuum 
permittivity; v 0 is incident light (laser) frequency; + and - signs refer to anti-Stokes and Stokes Raman 
scattering respectively; nk is the molecular vibrational frequency associated with the rotation-vibration 
transition from i to f, ga s the total degeneracy; £,is the energy of the initial state; k is the Boltzmann constant; 
T is the gas temperature; Q(T) is the total partition function of the molecule at temperature T; <£>(a , y, 6) is a 
function of the mean value of the derived polarizability tensor a 2 , the anisotropy y, and 0 is the observation 
angle that stands for the molecule dependent temperature-independent invariant, which is tabulated by Long 5 
for every branch of ro- vibration and pure rotation scatterings. Note that the factors (v+7 )b n 2 , and \b n 2 should 
be multiplied to Eq. (1) for ro-vibration Stokes and anti-Stokes scattering of diatomic molecule respectively, 
here v is the vibrational quantum number and b 2 is defined as h/87i 2 cnk where h is the Plank constant and c 
is the speed of light. The Raman spectral line locations, given by molecule vibrational frequency Vk, are 
calculated on the basis of energy conservation as 
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where Ef is the energy of the final state. The energy in a vibrational -rotational state (n, J) is expressed as the 
sum of the vibrational energy term, G(v) and rotational energy term, F(n, J) 
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According to the anharmonic oscillator model, G(v) of diatomic molecules may be given by 6 
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Here a) e , (*) e X e > an d ^eYe are molecular constants and tabulated by Hertzberg. 78 According to the non-rigid 
rotator model that takes the influence of centrifugal force due to vibration into account, the rotational term 
for the singlet state of diatomic molecules such as N 2 or H 2 may be given by 

F(y,J) = B V J(J + 1) - D v J 2 (J + l) 2 (5) 


where B n = B e - a e (n + 1/2), and D n = D e + P & (n + 1/2) where B t , D e , a, and /? e are molecular constants, 
which are tabulated by Hertzberg. 7 Thus for a Stokes Q branch, that is An = m-m = + 1 and AJ = //- /, = 0, 
Raman line frequency is given by 
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It should be noted that Q-branch of ro -vibrational Raman line is much stronger than the 0- and S-branch 
because the polarizability tensor greatly exceeds the anisotropy invariant. So we calculate only Q-branch of 
ro- vibrational Raman scattering for diatomic molecule in the present model. 

For S-branch (Stokes) of H 2 rotational Raman, the function ®(a 2 , y 2 , 0) is proportional to the 
anisotropy, y 2 and the following Placzek-Teller coefficients 5 
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Also, E(J) = F(J)/hc in Eq.(3) where F(J) = B e J(J+l) - D e J 2 (J+l) 2 , v k = E(Ji+2) - E(Ji) in Eq.(6), and 
Q(r) = Q r ot in Eq.(7) should be replaced. 


Thus, Raman line frequencies and cross-sections of ro-vibrational H 2 and pure rotational H 2 were 
calculated using above expressions. Note that we calculated the relative (but quantitative) intensity instead 
of absolute number by setting the function ® and molecular independent constants as a constant number so 
that the calculated Raman intensity can be fitted to experimental data of each molecule individually to 
determine temperature. 


The intensity I of the Raman lines follows from 9 


I = kngnll 0 (8) 

Here n is the molecular number density of the respective gas species in the measurement volume, Q is the 
solid angle observed, l is the length of the measurement volume, and I 0 is the laser intensity. Constant k is a 
setup constant and takes into account the geometry and wavelength-dependent losses. The differential Raman 

cross section — depends on molecule-specific parameters such as the isotropy and anisotropy of the derived 
polarizability tensor and also has a dependence on the geometry of the setup. For the configuration used here, 
a polarized excitation laser with polarization vertical to the scattering plane and observation of the scattered 
signal perpendicular to the laser propagation direction, the scattering cross section of the signal with the 
vertical polarization direction can be deduced by assumption of a Boltzmann distribution to be 
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For the horizontal polarization component one obtains 
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where C is a collection of several molecule-independent physical constants. For each vibration i, gi is the 
degeneracy; a{ and denote the mean values of the isotropy and the anisotropy, respectively, of the derived 
polarizability tensor. For linear molecules and in particular for two-atom molecules in general it is essential 
that 45(<2i’) 2 » 4(y /) 2 . For this reason the vertical polarization direction of the Raman signals of the vibration 
Q branch is more intense; i.e., the Raman effect essentially retains the polarization of the incoming laser 
beam. In contrast, fluorescence emission does not retain the laser polarization because of the relatively long 
time during which the molecule exists in the excited state before fluorescing. During this time the molecules 
rotate, eliminating any correlation between the polarization of absorbed and emitted photons. Thus, 
fluorescence signal contributions can be suppressed by use of a polarization sensitive detection scheme. 9 


A spectral profile of the Raman spectrum at atmospheric pressure is typically modeled using Gaussian 
profiles, neglecting collisional -broadening effects. 10 This assumption is even more valid at the low pressures 
of this investigation. However, both the laser and the spectrometer contribute to (broaden) the observed line 
shape of the Raman signal. The natural broadening laser and spectrometer convolve to form an overall 
Gaussian width of 11 

4v c (vJ) = (Avl + Avl aser + Avl vec r (! 1) 

Where, Av D is the molecular-dependent Doppler width, Av laser is the laser line width, 0.003 cm' 1 , and 
Av Spec is the spectral resolution of the spectrograph, 14.0 cm' 1 . For temperatures between 500 °C and 1000 
°C, Av Spec » Av d » Av Laser . The profile of Raman spectra at the temperature, T, was calculated by 
integration (sum) of all lines based on the assumption of the additive approximation. 

In order to quantitatively analyze the contribution of individual species in the overall Raman spectra, the 
relative scattering intensities of each of the different species had to be realized and properly implemented in 
the calculations. An empirical scaling factor for the scattering intensity of a rotational or vibrational Raman 
spectrum of species j relative to the H 2 O Raman scattering intensity can be defined as 12 : 

r Z5= a Ij /Nj(p,T) 
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Where Ij is the pixel intensity (CCD count per pixel) of a measured rotational or vibrational Raman spectrum 
of species j, n is the pixel number and a, b, c, or d indicate a corresponding pixel range for the summation of 
the spectrum of species j, and Nj(p, T) is the calculated chemical-equilibrium number density of species j for 
a given pressure and temperature. H 2 O was chosen as the normalizing species because its signal is free from 
the interference of the other major species in the plume. 

The data were reduced using a deconvolution routine written in Mathematica. While incorporation of 
dv 5pec was straightforward to shape the spectra for each rotational and vibrational transition, the masking of 
temperature dependence (captured in Av D ) precluded use of this model to estimate the temperature via line 
shape. Because of the relatively large spacing between the ro -vibrational H 2 Raman spectra (i.e. » Av Spec ), 
it was possible to estimate the (vibrational) temperature by fitting this spectra with peaks calculated via Eqns. 
3,5-7. The spectra lay between 665 nm and 690 nm, which is also largely free of overlap with Raman spectra 
from other major species. Within the error of measurement, the vibrational temperature will correspond to 
the kinetic/equilibrium temperature of the plume gas mixture. 


B. Sschlieren Diagnostcs 


Optical inhomogeneities refract light rays in proportion to their gradients of refractive index, n, in an 
x,y-plane. 13 Strong optical disturbances (such as shocks) will mask changes behind them. Depending on the 
degree of pressure matching, this may be true for Prandtl-Meyer expansion fans as well. However, the 
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boundary transition of the plume should yield data to approximate the velocity field if not the density field. 
The boundary of the plume can be characterized by its radius with respect to centerline and its local tangent 
angle, 0 = 8 + v e — v n , where 8 is the nozzle exit angle and v e and v n are the Prandtl-Meyer angles at the 
nozzle exit and after expansion, respectively. 14 8 and v n are constant for a given nozzle and O is taken from 
a Schlieren image. The local Mach number can then be determined from 15 



Unfortunately, the ratio of specific heats y, and the stagnation value of the density are unknown. They are 
approximated using the numerical modeling described below. 

The schlieren images were reduced using Matlab. Color images were first converted to greyscale for 
mathematical manipulation. A background (pre-firing) image was subtracted from each image taken during 
the thruster’s operation. Several images were then averaged and the contrast of the average was maximized. 
Plume angles, O, were determined by manual curve fitting of the final images. 

C. TDK- VIPER Modeling 

The Viscous Interaction Performance Evaluation Routine (VIPER) Version 3.6.2 (updated 2008) was 
used to model the expansion of the combustion gases from the exit of the catalyst bed to a slight distance 
downstream of the nozzle exit plane. The VIPER code is a modified version of the Two-Dimensional 
Kinetics (TDK) code with a Parabolized Navier Stokes (PNS) solver. 16 The modeling effort was undertaken 
to correlate the optically measured plume properties to pre -expansion (combustion chamber equivalent) 
values which could be used as inputs in the parallel effort to model the plume impingement. The model also 
provided a qualitative assessment tool for the concentration measurements. 

Initial temperature and pressure were taken to be experimentally measured values. The initial 
concentrations were adjusted to yield the observed concentrations in the near-field plume. No radial zoning 
was included as no variation in radial conditions could be tied to catalyst bed flow without post-test inspection 
of the catalyst bed and bed plate. Inert and slowly-reacting constituents were incorporated to simulate partial 
or off-nominal combustion. No liquid species were assumed. 

In principle, VIPER has the capability to extend the plume calculations downstream of the nozzle. 
However, specifying the appropriate mesh and viscous parameters proved impractical. Instead, a virtual 
nozzle corresponding to the plume boundary shown in the Schlieren images was used. This allowed 
calculation of the plume properties several nozzle diameters downstream. The impact of this approximation 
on centerline should be negligible, however, predictions near the edges of the plume will not include the 
transitional, viscous region of the real plume. 

III. Experimental Configuration 

This section describes the apparatus used to collect the Raman and Schlieren data. Brief descriptions 
are also provided for the thruster and test facility. 

A. GPIM Laboratory-model Thruster 

A lab model 22 N was provided by Aerojet-Rocketdyne for testing at NASA-GRC. The lab-model 
thruster incorporated a conical 8:1 nozzle with an exit pressure on the order of a few psia. This under- 
expanded condition in the RCL-chamber provided a target environment with identical plume constituents at 
lOx the pressure of the more flight-like expansion of a higher-fidelity thruster which has a pressure-matched 
nozzle. A more detailed description of the thruster is provided in an accompanying paper. 17 

Following a test series at Aerojet-Rocketdyne, Redmond, the thruster arrived at GRC with roughly half 
of its projected catalyst life remaining. Laser diagnostics did not begin until the thruster had completed 
additional acceptance testing at GRC, and, therefore, laser-based data collection began well into the mid-life 
of the thruster’s catalyst. 
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The thruster was typically operated in one of two modes during this test campaign. A pulse train of ten 
0.5 s pulses began and ended each test day as a thruster health check. The vast majority of the remainder of 
the segments were continuous 4 s firings with a few longer and a few shorter included to investigate specific 
trends. Chamber pressure and temperature measured downstream of the catalyst bed were recorded for each 
operating point as function of time at a rate of 4 Hz. This allowed spectroscopic and operational data to be 
correlated at several points during a firing segment. 

B. Altitude Simulation 

Plume measurements of an AF-M315E thruster were conducted at GRC’s Research Combustion Laboratory 
- Cell 11 (RCL-11). RCL-11 is a small rocket (<220-N), altitude (36.6 km) facility designed with optical 
access to incorporate the use of laser-based diagnostics. As shown in Fig. 1, the facility uses a six- foot long, 
three-foot diameter cylindrical vacuum tank. Thrust was measured via load cells. Vacuum is achieved by the 
use of a two-stage ejector system driven by motive air supplied from a central GRC facility. The thruster is 
fired horizontally into a water-cooled diffuser which provides an additional pumping effect for the vacuum 
tank. The water cooling system has a two stage dewatering subsystem to maintain the altitude condition of 
the vacuum tank throughout operation. RCL-11 has been used to test thrusters for up to 1 hour in duration, 
with the test duration limited by the thruster rather than the vacuum capability. 

The test environment in RCL-11 operated with a base pressure around 1720 Pa (0.25 psia) at quiescent 
conditions, and was seen to drop to 690 Pa (0.10 psia) during thruster firing. The chamber backpressure was 
primarily nitrogen gas, which was being used as a purge gas to cool the thruster injector between firings. The 
thruster was fired into a water-cooled diffuser with an inlet inner diameter of 8.26 cm (3.25 in) that was 
located approximately 6.35 cm (2.5 in) downstream of the thruster exit plane. Further details of the test cell 
are described in Deans. 18 

C. Raman Diagnostic 

As stated above, to circumvent issues associated with high-power laser delivery through vacuum 
chamber windows a high-frequency laser was coupled with an electron-multiplied ICCD (EMCCD) camera 
to mitigate the constraints imposed by high-power, single-shot lasers. A 10 kHz, frequency-doubled 
Nd:YAG solid state laser with 88 ns pulse widths yielded a time-averaged power of 200 W. The relatively 
low power-per-pulse enabled tight focusing of the beam through the vacuum chamber windows without risk 
of damage to the windows. The beam was delivered through a series of mirrors and collimating lenses from 
an adjacent laser room. Once at the test cell, the laser was directed downstream of the nozzle by a mirror 
and lens mounted on a translation stage which allowed remote positioning of the laser downstream of the 
nozzle. Figure 2 shows a schematic of the collection optics with respect to the test cell. 

Collection of the Raman signal was through a 6-fiber fiber optic cable. An imaging lens enabled each 
fiber in the cable to collect light from a separate collection volume focused 60 cm from the lens. The focal 
volumes were spaced 6.25 mm apart to simultaneously capture variations across the diameter of a pressure- 
matched nozzle. The under-expanded nozzle used in this investigation yielded a plume of smaller diameter 
and up to five of the interrogation points lay outside of the plume depending on axial location. This allowed 
simultaneous measurement of the plume and the ambient gas, shot-to-shot calibration of the laser power, 
detector drift, and laser-detector alignment. Measurement of signals outside of the plume also provided a 
path to eliminate facility effects such as window fluorescence. A three stage positioning system remotely 
translated the collection optics and allowed optimization of the laser-collection optics alignment at 
atmospheric pressure. Negligible hysteresis was observed in the stages which allowed rapid and repeatable 
interrogation of locations downstream of the nozzle exit plane. 

The multi-fiber cable delivered the Raman scatter to a back-illuminated EMCCD camera via an imaging 
spectrograph. The PI-MAX4:512EMB used in this investigation fiberoptically coupled the EMCCD to an 
image intensifier that enabled six times higher light throughput than lens -coupled configurations. The camera 
potentially allows lOOOx greater sensitivity than standard ICCD cameras. 16 In addition the intensifier can be 
gated to the pulse width of the laser at the laser frequency while the EMCCD accumulates tens to over 10000 
laser-shots to form a single image. In preparation for testing on the GPIM thruster, Kojima demonstrated 
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this technique’s ability to capture time -resolved characteristics of an atmospheric hydrogen flame. 19 The 
imaging spectrograph allowed each of the fiber’s signals to be independently imaged on the camera and thus 
simultaneously collected through independent accumulations. 

A sharp-edge filter suppressed the resonant (Rayleigh) scatter from the plume and all anti-Stokes Raman 
spectra. Raman spectra from just above 532 nm to 700 nm was collected on the camera. This range was 
sufficient to include rotational and vibrational spectra. Spectral resolution was sacrificed in order to 
maximize the throughput (signal) of the collection optics. This introduced significant instrument broadening 
of the Raman signals on the order of several nm. To the extent possible, this broadening was removed in the 
deconvolution as described above, but significant uncertainty remains due to the overlap of N 02 ,N 20 , and 
CO 2 vibrational and H 2 rotational spectra in the 550 nm to 585 nm range. 

D. Schlieren Diagnostic 

A coincident configuration was selected for the Schlieren diagnostic to facilitate rapid, remote transition 
from the schlieren to the Raman diagnostic. This technique utilizes a single primary mirror instead of two as 
in the common z-type schlieren optics configuration. A schematic of the coincident Schlieren diagnostic 
configuration is shown in Figure 3. The coincident configuration offers several advantages including a double 
pass through the target area which, in principle, doubles the sensitivity of the measurement and reduction in 
experimental footprint with long-focal length mirrors. 

Background-oriented Schlieren has been demonstrated to provide the best quantitative analysis of a 
similar flow field. 13, 20 However, this method was not incorporated because it would have precluded viewing 
of the thruster by facility monitoring cameras and because the background may have fluoresced during 
Raman interrogation. 

The large turning mirror introduced distortion, but atmospheric and subsequent vacuum data suggested 
that this effect could be neglected in the determination of the plume shape. The ribbon-filament tungsten 
lamp provided a wide-band, variable intensity light source. Optics incorporated in the configuration 
eliminated the ultraviolet and far-infrared light. A 50-cm diameter, 1 -meter focal length spherical mirror was 
used as the primary mirror. A larger diameter mirror would have been preferable, but it would have interfered 
with both the Raman diagnostic and optical access in general to the thruster. As shown in the schematic, 
only the primary mirror was located in the vacuum chamber. A low distortion, F4.5-F22C zoom lens with 
a relatively large depth of field was used to image a region slightly larger than the primary mirror. 
Optimization of the coupling lens yielded a magnification of roughly 3:1. Baffles were employed to prevent 
stray light from the tungsten lamp from entering the camera. 

IV. Results and Discussion 

Results of the Raman and Schlieren diagnostics are discussed below followed by a correlation of those 
measurements with the VIPER modeling. The laser had a long beam waist which allowed for collection of 
data in the vertical direction without correction for laser intensity. This was verified repeatedly during 
alignment checks which are described below. 

The thruster was operated in a manner that facilitated plume interrogation and not in ones that 
corresponded to either design operation: pulse trains of multiple, short duration pulses or steady-state 
operation. Typical firings consisted of 4-s pulses which approached steady conditions in pressure and flow 
near the 4-s point and thus pointed to concentrations that might be expected in truly steady operation. A pair 
of 8-s firings were included in the test matrix to verify that the concentrations seen at 4 s were representative 
of steady operation. Pulse trains of ten 0.5 -s pulses were undertaken as daily health checks of the thruster. 
As the test campaign progressed, these provided an opportunity to interrogate the plumes associated with this 
mode of operation as well. 
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A. Typical Raman Spectra 


At the beginning and end of each test sequence, Raman data were collected at atmospheric pressure to 
calibrate the spectra given changes in laser power, window and/or mirror degradation, and spectrograph drift. 
These calibration checks also allowed for optimization of the Raman signal collection via alignment of the 
collection volumes with the laser. Data taken post-test at roughly atmospheric pressure had very weak O 2 
and H 2 O signals as the rise in pressure was mostly due to the N 2 facility purge. The two “high-pressure” 
conditions facilitated calibration and validation of the deconvolution model. Raman spectra were taken at 
simulated altitude conditions without thruster operation. Similar, slightly weaker, spectra were obtained 
outside of the plume during thruster operation. These vacuum, no plume data enabled calibration and ensured 
alignment of the laser and collection volumes. They also facilitated proper adjustment of camera integration 
and gain settings. 

Figure 4 shows a Raman spectrum taken during thruster firing at catalyst mid-life with p c = 1.55 MPa 
(225 psia) and Tc = 1480 K. As indicated in the figure, major plume species associated with equilibrium 
combustion H 2 O, N 2 , H 2 , and CO 2 are present. The spectrum shown in Fig. 4 resulted from the accumulation 
of 1000 individual spectra (0.1 s total accumulation time). A fit also is shown which is used to deconvolve 
the various individual spectra. From Eqns. 9 and 10, if the laser’s polarization is shifted 90°, the resulting 
“Raman” spectrum should not contain vibrational -Raman spectra. Figure 5 shows the spectrum resulting 
from the subtraction of the polarization shifted signal from that in Fig. 4. This spectrum should contain only 
vibrational-Raman spectra and therefore is better- suited to identifying the major species in the plume. 

Repeated collection of atmospheric Raman signals eliminated uncertainty in the ICCD pixel -to- 
wavelength correlation. It also removed uncertainty associated with large, day-to-day variations in laser 
quality (which were small). Impacts due to shot-to-shot variations in laser power were mitigated by 
accumulating large numbers (100’s) of individual shots for each data point and, when possible, averaging 
several data points to obtain a single spectra. No significant variation in laser power was observed over the 
course of a day. 

However, uncertainties remain in the amplitude of individual peaks in the spectra, in the species 
associated with peaks in regions of Raman signal overlap, and in the cross sections associated with many of 
the vibrational Raman spectra necessary for estimating concentrations. Vibrational N 2 , H 2 O, and H 2 Raman 
spectra present in regions relatively free of other spectra. In addition, the signal to noise ratio for these 
transitions is greater than 5 for most cases. Therefore, the uncertainty in estimating peak height associated 
with these species is less than 10% for N 2 and H 2 O and 20% for H 2 . As seen in Table 1, the spectra of CO 2 , 
NO 2 , N 2 O, and N 2 O 2 lay in the same region and, in part due to instrument broadening, convolve to yield a 
series of weak peaks with significant uncertainty in the individual contributions of each of these species 
overall spectra. While tempting to give preference to CO 2 as the only equilibrium combustion product in this 
group, the presence of hydrocarbons in the spectra suggest it is unlikely to be present in concentrations 
associated with equilibrium. It is impossible a-priori to say which species is more prevalent and non-unique 
deconvolutions result which match the spectra with 30-40% uncertainty in each of these species’ 
concentrations. Because the N 2 O cross sections are so large relative to the others, the relative similarity in 
overall peak strengths suggest its concentrations are always relatively small, so the impact of a large 
uncertainty in its concentration does not correspond to a similar uncertainty in the concentrations of other 
species. 

Temperature at the point of interrogation was estimated both by fitting the N 2 rotational peak and the H 2 
ro- vibrational spectra. The N 2 -based estimate was limited by the instrument broadening and by the presence 
of the large non-Raman signal at 545 nm. The latter was removed via a Gaussian fit to its spectra. 
Uncertainty in the temperature measurement increased during the test campaign from roughly 10% due only 
to instrument broadening to 30% at the end of testing due to the significant convolution of the N 2 and 545 
nm spectra. For those conditions in which signal to noise of the H 2 was greater than 3 (most of the data after 
3 s during a thruster firing), fitting of the H 2 ro- vibrational yielded temperature estimates comparable to the 
N 2 fits. While this was reassuring, it did not reduce the overall uncertainty in temperature because the 
sensitivity of the ro- vibrational spectra could not resolve temperature changes below roughly 15%. At the 
end of testing, when the N 2 -based temperature measurement had an uncertainty greater than this, the H 2 
spectrum was reduced in magnitude because of the increase in hydrocarbon capture of the hydrogen. 
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B. Trends in Concentrations 


The Raman diagnostic was intended to measure the species concentrations and temperature at a couple 
of steady-state, representative operating conditions and to capture trends in the specie concentrations with 
chamber pressure, time on catalyst, and as a function of position in the plume. Each of these is discussed 
below. 

Figure 6 compares the measured concentrations for chamber pressures of 170, 200 and 250 psia taken 
3 mm downstream of the nozzle exit plane at catalyst mid-life. Raman data were accumulated for 0.2 s and 
three data sets were acquired and averaged over the last second of the 4-s firings at these conditions. 
However, except for the increase in Xh 2 o, all changes in concentrations are within the uncertainty of the 
measurements and suggest that the pressure has little impact on plume constituents over the ranges tested. 

Figure 7 shows the variation in specie concentrations during an 8-s firing during catalyst mid-life with 
Pc = 230 psia for the final 5 seconds. As alluded to in the above discussion, there was no discernable Raman 
spectrum in first 1 .5 s of operation. Raman data were accumulated for 0.2 s and four data sets were acquired 
and averaged over each second of the last 7 s of the 8-s firing. (Data points in Fig. 7 represent the centers 
of these averages.) There are no significant changes in any trend after 4 s which validates the 4-s firing as a 
minimum firing time to reach steady conditions. Between 1.5 and 3.5 s, there are significant changes in the 
mole fractions of the constituents. 

Figure 9 shows spectra taken simultaneously at 6 radial locations 6.25 mm apart 4-mm downstream of 
the exit plane for Pc=250 psia. Data were accumulated for 1 s on the ICCD and the spectra falling roughly 
between 2.75 to 3.75 s were used to approximate the steady-state concentrations. Fong integration times 
were required because larger gains associated with fewer accumulations resulted in cross-talk between the 
fiber signals. Rotational and fluorescence data have not been removed since three of the interrogation points 
lie outside of the plume and are dominated by rotational-N 2 spectra. Figure 10 compares the calculated mole 
fractions for these locations. Variations in the concentrations are within the uncertainty of the calculations. 
Figure 1 1 shows a slight variation in temperature with radius. Note the temperature outside of the plume 
suggests some mixing, but the instrument broadening introduces larger error at lower temperatures. 

Figure 12 shows the variation in concentrations as a function of axial position along centerline. There 
is very little variation as would be expected in the high-speed (frozen) flow. Temperature variation with 
axial position is given in Fig. 13. Temperatures calculated using rotational-N 2 spectra and ro-vibrational-H 2 
spectra are compared in the figure. As the signal decreases away from the exit plane, the H 2 -based 
temperature has very large uncertainty. However, at 1 and 3 mm, the temperatures agree within the 
uncertainty. 

To complete the discussion of concentration trends in the catalyst mid-life, Figure 14 shows the typical 
variation in concentrations during a pulse train of 10 0.5 -s pulses. Data were accumulated for 0.1 s and the 
data falling roughly between 0.35 and 0.45 s in a pulse were used to characterize the pulse. Trends seen 
during the multi-second firings above were also observed in both the pulse train and individual pulses. The 
first pulse was dominated by a non-Raman signal. The next three pulses have large non-Raman signals, but 
Raman signals are present and relative concentrations of the major species can be approximated. The final 
six pulses have clean Raman signals and yield quantitative mole fractions. Note that there is a monotonic 
trend in the CO 2 with pulse number. It corresponds to the trend in Tc as shown in Fig. 15 and to some extent 
the trends in Fig. 8. Each pulse had a similar non-Raman start and finish as the longer 4-s firings, though 
both were smaller for the shorter firings. 

Pulse trains at catalyst end-of-life showed similar trends to those at catalyst mid-life with the exception 
that the fluorescent signals remained throughout the pulses and occasionally a mid-train pulse exhibited 
behavior similar to the first pulse. Figure 16 shows the trends in concentrations as measured by 0.1 s 
accumulations between 0.35 and 0.45 s of 0.5 s pulses. 

The Raman diagnostic proved successful in measuring both specie concentrations and temperatures in 
the plume of the 22-N laboratory-model thruster. The Raman diagnostic technique of using a high-frequency, 
relatively low power laser and accumulating many individual spectra on an EM -ICCD was demonstrated to 
be capable of characterizing the plume of the GPIM lab -model thruster. High signal-to-noise data were 
collected of the background gas in the test chamber which was roughly at the pressure expected in the plume 
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of the EM thruster. This suggests that the Raman technique has the sensitivity to measure the specie 
concentrations in the plume even at pressures one-tenth that seen in the present investigation. 

C. Schlieren Imaging 

Schlieren imaging was implemented at the end of the test campaign. For assessment of plume structure, 
the Schlieren images were integrated for 50 ms and taken in a continuous stream from pre-firing to post- 
firing. 50 ms integration provided a compromise between intensity and time resolution. Longer tended to 
blur the transition region due to averaging of oscillations in the plume boundary (likely corresponding to 
oscillations in Pc) and shorter exposures yielded weaker overall signals. 

Typical raw and reduced Schlieren images are shown in Figure 17. The image is of 3.5 s into a 4-s firing 
with Pc=200 psia. Note that there are clear boundaries associated with the viscous layer extending roughly 
5 mm downstream in the image. A core flow is evident as are regions of large gradient near the exit of the 
nozzle. (The bow- shock- like structure on the right hand side of the image is likely a false image resulting 
from contamination on one of the mirrors as the structure is present in images without flow. In subtracting 
the background image, it presents as a white structure.) 

A series of Schlieren images in Fig. 1 8 shows the development of the plume over a 4-s firing with Pc=200 
psia. Note that between 0 and 1.3 s there is a strong plume boundary suggesting reacted, supersonic flow 
and also a strong core flow. This supports the hypothesis that the large fluorescent signal at the beginning of 
firing results from gaseous, incomplete combustion including a large amount of soot-like effluent. As the 
firing progresses, the plume boundary is seen to weaken which is largely due to the large pressure oscillations 
occurring at this point in the test campaign. At the end of the firing, the flow resembles the pre-firing flow 
and there is no structure that would support a large scattering from uncombusted propellant. This supports 
the hypothesis that the end-of-firing fluorescence results from a diffuse flow back into the test chamber from 
the diffuser. 

Figure 19 shows velocity profiles calculated from Schlieren images for operation at Pc = 110 and 200 
psia. There are large uncertainties in the measured values resulting from the images being slightly out of 
focus. The trends and magnitudes are, however, consistent with those calculated using the VIPER code. 

D. Correlation of Data to VIPER Modeling 

Flow characteristics upstream of the nozzle exit plane calculated using the VIPER code to match the 
concentrations and temperatures measured via the Raman technique in the plume. The chamber pressure and 
temperature were set using experimentally measured values as well. The relative concentrations of the major 
species are more consistent with equilibrium flow than frozen flow, their values being reduced by the 
additional carbon, methane, and nitrogen-oxides present. However, the VIPER-predicted plume 
temperatures were consistently about 300 K less than those measured. Plume temperatures calculated using 
the HAP code and identical chamber conditions, 17 however, agree with those measured. 

Unfortunately, the large uncertainty in mass flow rate through the thrust chamber precluded comparisons 
of predicted and measured performance. However, trends in the predicted and measured performance were 
consistent. In particular, inclusion of soot and larger hydrocarbons in the gas constituents resulted in 
performance degradation similar to that seen between mid-life and end-of-life operation. 

No effort was made to capture the large fluorescent signals during the start of each firing. However, 
soot and small hydrocarbons were incorporated upstream of the throat of the nozzle with some success. This 
suggests that VIPER can be used to model the performance of AF-M315E thrusters given reliable predictions 
for the plume species as a function of chamber pressure. This investigation provides an advance in that 
direction. However, higher-resolution Raman data are required at or before catalyst mid-life to support 
mission planning. 


10 



IV. Future Work 


The diagnostic effort reported in this publication was intended to be a preliminary assessment of the 
diagnostic techniques for subsequent implementation on a higher fidelity AF-M315E fueled thruster. While 
the second round of testing was cancelled, several lessons learned in this first test campaign boded well for 
a more robust characterization of the EM thruster. 

A. Improvement in Raman Characterization 

Reduction in the instrument broadening which led to large uncertainties in both the amplitudes and 
associated species of Raman signals in the 550-580 nm region can likely be achieved by replacing the grating 
in the spectrograph. Because of the large signal-to-noise in the background gas, it is projected that the loss 
in signal with the finer grating can be tolerated. This improvement will allow discrimination between the 
nitrogen-oxides (typical of HAN combustion/decomposition) and hydrocarbons (typical of HEHN) both 
when combined in AF-M315E and when tested separately. 

B. Improvement in Schlieren Characterization 

Introduction of a higher-fidelity knife edge resulted in a significant gain in signal to noise ratio of the 
Schlieren diagnostic. In addition, a simple optimization of the alignment removed distortion evident in Figs. 
17 and 18. A larger primary mirror was ordered for use in the 2 nd series of tests. 

Longer integration times were projected for interrogation away from the end-of-life condition. This too 
will greatly enhance the signal to noise ratio and will likely enable the technique to be applied to the higher 
expansion ratio nozzle. However, as the nozzle would likely have been pressure matched, the transition from 
plume to background would have been much more subtle. 

V. Summary and Conclusion 

The novel approach of accumulating a large number of individual Raman spectra generated by a high- 
frequency laser was demonstrated. This was enabled by the new capabilities of the electron -multiplied ICCD 
camera. Raman characterization of the plume of a 22 N laboratory model thruster showed concentrations of 
major species consistent with predictions based on equilibrium conditions with the exception of added 
hydrocarbon species. Strong N 2 rotational and vibrational Raman signals were obtained even at pressures 
below 0.1 psia. Significant improvement in the spectral resolution is believed possible through the simple 
changing of the discriminating grating in the spectrograph without prohibitive losses in signal strength. 

Schlieren imaging provided the first structural characterization of an AF-M315E plume. The images 
indicated plume velocities consistent with VIPER model predictions. However, the resolution of the images 
was not sufficient to enable quantitative characterization of the 2-dimentional density gradient in the plume. 
Post-test replacement of the knife-edge yielded much sharper imaging capability which should enable this 
quantitative characterization in future tests. Qualitative analysis of the plume structure showed the schlieren 
technique’s ability to capture plume oscillations on the timescale of 50 ms. 

The VIPER code was successfully modified to predict the kinetic expansion of the AF-315E 
monopropellant through the thrust chamber and downstream of the nozzle. The latter was accomplished by 
matching the expansion to the boundary measured via schlieren. Hydrocarbons and nitrogen oxides were 
included in the non-equilibrium composition and allowed to react. Very little change in concentrations during 
the expansion were observed. 

The test campaign was successful in generating sufficient data to anchor the plume impingement model 
presented in a companion paper. Through straight-forward modifications, the techniques appear quite 
capable of characterizing the plumes of both higher-thrust and higher expansion nozzles of thrusters using 
green-propellants such as AF-M315E or other HAN and HEHN based propellants. 
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Table 1 Parameters of Raman Spectra for Molecules of Interest. 6 ’ 8 ’ 9 ’ n Error! Bookmark not defined. 
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Figure 1 Schematic of the RCL-1 1 test chamber. 



Figure 2 Schematic of laser delivery to the test facility. 




Figure 3 Schematic of Schlieren optical configurations. 13 The single-mirror coincident configuration was 
used in this investigation. 
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Figure 4 Typical mid-life Raman data and fit. 
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Figure 5 Vibrational-Raman signals remaining when the non-polarized signal (rotational-Raman and LIF) is 
subtracted from the total collected signal (without half-wave plate). 
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Figure 6 Trends as a function of chamber pressure, Pc, at catalyst mid-life. 
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Figure 7 Variation of major species during an 8 s pulse at catalyst mid-life. 
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Figure 8 Trends in plume temperature and chamber temperature during an 8 s pulse at catalyst mid-life. 
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Figure 9 Spectra collected simultaneously at 6 different radial locations 5 mm downstream of the nozzle exit 
plane. 
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Figure 10 Concentrations at 6 radial locations based on the spectra in Fig. 9. 
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Figure 1 1 Radial variation in temperature measured 5 mm downstream of the exit plane. 



Figure 12 Variation in concentrations as a function of distance downstream of the nozzle exit plane. 
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Figure 13 Variation in temperature as a function of distance downstream of the nozzle exit plane. 
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Figure 14 Variation in concentrations during a pulse train at catalyst mid-life. 



Figure 15 Variation in chamber temperature during a pulse train at catalyst mid life. 
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Figure 16 Variation in concentrations during a pulse train near catalyst end-of-life. 



Raw schlieren image 



Figure 17 Typical Schlieren raw and enhanced images. 
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Figure 18 Series of schlieren images showing variation in plume structure over a 4-s firing near catalyst end 
of life. Pc = 200 psia. 
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Figure 19 Velocity as a function of distance from nozzle exit plane calculated from schlieren images for 
chamber pressures of 1 10 and 200 psia. Calculated values using the VIPER code are also shown. 
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